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We derive a new class of exact solutions of Einstein's equations providing a physically plausible 
hydrodynamical description of cosmological matter in the radiative era, between nucleosynthesis and 
decoupling. The solutions are characterized by the Lemaitre-Tolman-Bondi metric with a viscous 
fluid source, subjected to the following conditions: (a) the equilibrium state variables satisfy the 
equation of state of a mixture of an ultra-relativistic and a non-relativistic ideal gases, where the 
internal energy of the latter has been neglected, (b) the particle numbers of the mixture components 
are independently conserved, (c) the viscous stress is consistent with the transport equation and 
entropy balance law of Extended Irreversible Thermodynamics, with the coefficient of shear viscosity 
provided by Kinetic Theory. The fulfilment of (a), (b) and (c) restricts initial conditions in terms 
of an initial value function, A' s \ which in the limit of small density contrasts becomes the average 
of spatial gradients of the fluctuations of photon entropy per baryon in the initial hypersurfacc. 
For A^ s ' 7^ and choosing the phenomenological coefficients of the "radiative gas" model, we have 
an interactive photon-baryon mixture under local thermal equilibrium, with radiation dominance 
' and temperatures characteristic of the radiative era (I0 6 K > T > 10 3 K). Constraints on the 

observed anisotropy of the microwave cosmic radiation and the condition that decoupling occurs 
at T = T D w 4 x I0 3 K yield an estimated value: jA^ s) | » I(T 8 which can be associated with a 
bound on promordial entropy fluctuations. The Jeans mass at decoupling is of the same order of 
magnitude as that of baryon dominated perturbation models (« 1O 16 M0). 
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I. INTRODUCTION 
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The radiative era of cosmic evolution comprises the period from the end of primeval nucleosynthesis to the decoupling 
(j \ of matter and radiation (see refs 0] to Ell). A gross description of cosmological matter sources in this period is given 
CJ^' by an interactive mixture of ideal relativistic and non relativistic gases ("radiation" and "matter") in local thermal 
^ ( | equilibrium (LTE). 

The standard approach to this type of matter source is either a FLRW spacetime whith equilibrium Kinetic Theory 
distributions ||, ||, 0, gauge invariant perturbations on a FLRW background Q, ||, p0[ , and various 
types of hydrodynamical models, fu^j , pifl , (m) which, in general, fail to incorporate a physically plausible description 
of the interaction between matter and radiation. Even if we argue that the universe is "almost FLRW" or "almost 
in thermal equilibrium", the small deviations from equilibrium are extremely important, [Q, ||, 

to account for most interesting phenomena of cosmic evolution: nucleosynthesis, structure formation, abundance of 
relic gases, etc. Models with perfect fluid sources, whether hydrodynamic, |K| , [ jl3| , or based on Kinetic Theory, [ jl5| , 
necessarily assume a quasi-static adiabatic and reversible evolution and thus, fail to incorporate into the resulting 
picture even small deviations from equilibrium. 

Dissipative sources have been incorporated numerically within a purely FLRW geometry [[l6f or following a pertur- 
bative approach jl?]]. However, the literature still lacks an alternative hydrodynamical treatment, based on inhomo- 
geneous exact solutions of Einstein's equations with dissipative sources and fully complying with the thermodynamics 
of a radiative gas within a transient regime. Ideally, such exact models should include all dissipative agents (heat 
flux, bulk and shear viscosity) and should be consistent with the theoretical framework of Extended Irreversible 
Thermodynamics (see refs JlT[ | to p?|), thus satisfying suitable transport equations complying with causality, with 
phenomenological coefficients given by Kinetic Theory for this type of source. Since this general treatment would be 
mathematically untractable, we aim at the best possible approach based on exact solutions of Einstein's equations. 
Therefore, we have made the following simplifying assumptions: (a) the matter source is a fluid with shear viscosity 
but without heat conduction nor bulk viscosity, (b) the equilibrium state variables satisfy the equation of state of a 
mixture of relativistic and non-relativistic ideal gases, where the internal energy and pressure of the latter have been 
neglected, (c) the particle numbers of each mixture components is independently conserved, (d) we exclude dark mat- 
ter and/or exotic particles and assume instead a tight coupling between photons (radiation) and baryons and electrons 
(matter), hence there is a common temperature for the mixture (LTE), while the microscopical interaction models 
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are the various processes of radiative transfer, [Q, ||, {I, p3[ , [ p5| , (26), |27|: Thomson scattering, brehmstralhmg, 
free-free absortion, etc. Although this type of interactions involve mostly photons and electrons, the dynamics of the 
matter component is governed by the baryons since the latter provide most of the rest mass content of non-relativistic 
matter (without dark matter). 

Restrictions (b) and (c) are easy to justify: since the ratio of photons to baryons is such a large number 10 9 ), 
we can truly ignore the pressure of non-relativistic matter. Also, after nucleosynthesis, in the temperature range 
T < 10 6 , matter creation and anhilitation processes balance each other and effectively cease to be dynamically 
important ||, ||. On the other hand, the lack of heat conduction (restriction (a)) is more difficult to justify. It 
can be associated with an adiabatic (zero heat flux) but still irreversible evolution (nonzero viscosity), and can be 
a reasonable approximation on specific conditions. For example, for a radiative gas at higher temperatures shear 
viscosity dominates over heat conduction, but the latter becomes significant as the mixture cools [g]. The lack of bulk 
viscosity is a better approximation: it is negligible for a radiative gas in the temperature range 10 3 K < T < 10 6 K that 
we are interested, PL M L P l L 23 1 , p5[ , becoming important for higher temperatures (the mid relativistic regime 
where k B T rj mc 2 fll7|], p8[, |]29|] ) . However, we accept that ignoring these dissipative fluxes weakens the scope and 
validity of the models, but we argue that this is compensated by the simplification of the field equations, leading to 
exact forms for the equilibrium state variables and shear viscosity that still satisfy (under the restrictions mentioned) 
thermodynamically consistent relations. 

The models we present are based on the spherically symmetric Lemaitre-Tolman-Bondi metrics, usually associated 
with dust sources f3Cf| , pl|| . However, this metric is compatible with a comoving fluid source with zero heat flux but 
with anisotropic stresses, which we describe as shear viscosity. Obviously, the lack of heat flux and 4-acceleration 
necessarily implies a very special shear viscous tensor whose divergence exactly balances the nonzero spatial gradient 
of the equilibrium pressure. Considering this metric and this source, we impose on the equilibrium state variables 
the equation of state for a mixture of ideal gases (under the restriction (b)). The field equations can be solved up 
to a quadrature, without having to make any assumption on the form of the shear viscous pressure. The latter, as 
well as all equilibrium state variables can be determined from the solution of the quadrature, up to two initial value 
functions that can be identified with the initial energy densities of the matter and radiation components. We consider 
only the case that would be equivalent to spacelike sections of zero curvature. A generalization of this class of exact 
solutions to the more general Szekeres-Szafron metrics admitting no isometries has been published recently p2| , while 
the study of a non-relativistic ideal gas is considered in J33| . 

Once the field equations have been integrated, we define a set of initial value functions that gauge the deviation 
from homogeneity of the average of initial density contrasts. The terms involving various gradients of metric functions 
can be given in terms of these gauges, so that in the limit when the latter vanish a FLRW spacetime can always be 
obtained as the homogeneous (and reversible) subcase. In section VII we derive the conditions that the models must 
satisfy in order to be consistent with the theoretical framework of Irreversible Extended Thermodynamics, in the case 
where shear viscosity is the only dissipative agent and the coefficient of shear viscosity is that given by Kinetic Theory 
for the radiative gas ||, 23 1, pq , p7j , p9| . This leads to an entropy balance law and a suitable transport 
equation for shear viscosity that is satisfied for a specific functional form of the relaxation time. Conditions are given 
so that the latter quantity behaves as a relaxation parameter for an interactive cosmological mixture of matter and 
radiation. These conditions of thermodynamical consistency are then explicitly tested on the models, leading to a set 
of restrictions on the initial conditions (the latter given in terms of the gauges of initial density constrasts). The most 
relevant result is that thermodynamical consistency constrains an initial value adimensional function, , which in 
the limit of small density contrasts is approximately the average gradient of the photon entropy per baryon along the 
initial hypersurfacc t = ti. An analogy is provided with the theory of perturbations on a FLRW background, whereby 
A> s ' = is formally analogous to the definition of initialy adiabatic perturbations in the sincronous gauge Q], [[|, 
||, ||lo| . The constraints on the observed anisotropy of the microwave cosmic background, as well as the condition 
that decopling occurs at T = T D « 4 x 10 3 K, leads to the estimated value |A| S '| 10~ 8 . Since initial conditions of 
the radiative era should be traced to previous periods of cosmic evolution, this constraint can be related to maximal 
bounds on entropy fluctuations in primordial perturbations. Finally, we compute the Jeans mass associated with the 
thermodynamically consistent models, leading to a value similar to that obtained for baryon dominated perturbation 
models: M, w 10 16 M n . 



II. INTERACTING MIXTURE OF RADIATION AND NON-RELATIVISTIC MATTER 

A radiation-matter mixture can be described by a mixture of two ideal gases: one an ultra-relativistic gas of massless 
particles, the other a non-relativistic ideal monatomic gas with m being the mass of the particles. This is characterized 
by the total matter energy, p and pressure p 
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p = mc 2 n {m) + ^n {m) k B T {m) + 3n (r) fc B T (r) , (la) 

p = n (m h B T^ + n {r ^k B T {r \ (16) 

where k B is Boltzmann's constant and n,T are particle number densities and temperatures of the two components, 
distiguished by the superindices (m) ("matter") and (r) ("radiation"). If there is local thermal equilibrium (LTE) 
between the components, the latter interact and evolve with the same temperature: — T^ m%> = T. If the 

components are decoupled, each gas evolves with a different temperature. 

Assuming LTE, if <C n^~\ but the ratio mc 2 /k B T is not negligible, then equations (1) can be approximated 

by 

p a mcV" 1 ' + 3n (r) fc B T, (2a) 

P«« lr) W (26) 

an equation of state describing a radiation dominated mixture in which the presence of non-relativistic matter is 
dynamically important. If we assume non-relativistic matter to be made up of baryons (with m being a protonic 
mass) and since the ratio of baryons to photons n^/n^ w 10~ 9 is a small number, the equation of state (2) is 
a reasonable approximation in the temperature range 10 3 < T < 10 6 K, characteristic of the "radiative era" from 
the end of nucleosynthesis to the transition between radiation to matter dominance, including the recombination 
and decoupling eras. At such temperatures, it is also safe to assume, jy, ||, j|, ||, Q that electrons and photons 
interact mostly through Thomson scattering but creation and annihilation processes (bremsstrahlung and free-free 
absorption) roughly compensated one another so that particle number densities of the components of the mixture 
satisfy independent conservation laws. Once the decoupling of the matter-radiation mixture takes place at about 
T w 4 x 10 3 K, the assumption of LTE is no longer valid and interaction between components ceases. Equation of 
state (1) can also be approximated by a form similar (2) with the internal energy of radiation taking approximately the 
Stefan-Boltzmann law: p^ — a^T 4 , where as denotes the radiation constant. However, out of thermal equilibrium 
the Steffan-Boltzmann law is incompatible with the ideal gas equation of state. 

Having in mind the conditions justifying (2), we will describe a matter-radiation mixture evolving along adiabatic 
but irreversible processes by the fluid tensor 

T ab = pu a u b + ph ab + n Qb , (3) 

h ab = c- 2 u a u b + g ab , u a Ii ab = 0, n a a = 0, 

where: p,p satisfy (2), u a is the 4- velocity shared by radiation and matter, IP 6 is the shear viscous pressure tensor 
(a symmetric traceless tensor) which arises because of the matter-radiation interaction, and particle number densities 
satisfy the conservation laws 

(n< m V) ;a = 0, {n^u a ), a = 0. (4) 

As mentioned previously, bulk viscosity is negligible within the temperature range we are interested in pah Q , |25| , 
p6| , |27|] , |2§f| , while even if neglection of heat conduction can be justified for relativistic temperatures||] , it does 
weaken the scope of the models. However, this restriction is compensated by the obtention of exact solutions that are 
still thermodynamically consistent. 



III. THE LEMAiTRE-TOLMAN-BONDI METRICS 

Consider (3) as the source of the "Lemaitre-Tolman-Bondi" (LTB) metric ansatz, usually associated with spherically 
symmetric Lemaitre-Tolman-Bondi dust solutions |50|, |3l]] 
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ds 2 = -c 2 dt 2 + dr 2 + Y 2 (d0 2 + sm 2 (6)d(j) 2 ) , (5) 

1 — F 

where Y = Y(t, r), F — F(r), and a prime denotes partial derivative with respect r. Just as in the LTB dust solutions, 
we assume the coordinates in (5) to be comoving and the 4- velocity of the fluid source to be u a = cS a t, a geodesic 
vector field, since u a = u a -bU b = 0. Other kinematic invariants associated with (5) are the scalar expansion O = u a ;a , 
and the shear tensor a a b = u^ a . b ^ — (0/3)/i a h, given in the coordinates of (5) by 

V' 2Y 

6 = — + — , (6a) 
Y> Y 



a a b = diag [0, -2a, a, a], a = - ^- - — J , (66) 

while the most general form of II a b for the metric (5) is given by 

IF 6 = diag [0, -2P, P, P] , (6c) 

where Y = u a Y^ a — Yj and P = P(t, r) is an arbitrary function. Notice that a comoving and non-accelerating 
4- velocity does not imply p' — 0, as in the perfect fluid case (G e q — G r r — 0). As revealed by the momentum balance 
law: h ca T ab = 0, applied to the viscous fluid source (3), we have 

h a b {p,b +H bc;d h cd )=0 (p - 2P) ' + 6Py = (7a) 

showing how the divergence of the shear viscous tensor exactly balances the nonzero pressure gradient. The energy 
balance, u a T ab -b = 0, is given by 

p+ (p + p)Q + a ab U ab = => p+^Qp + 6<TP = (76) 

illustrating how the term a a bH ab = 6trP can be understood as an interaction term responsible for local energy exchange 
between matter and radiation. 

Integration of the conservation laws (4) for (5) yields 

= f$V n« = ( $ V , (8) 

1 \Y j Y'/Y' 1 \Y ) Y'/Y ' v ; 

where n[ m \ n| r ' depend only on r and are the particle number densities of non-relativistic matter and radiation, 
evaluated along a suitable initial hypersurface labeled by t = U. The subindex affixed to any quantity, as Y"i, will 
denote henceforth initial value functions (functions of t,r evaluated along t — ti). It is important to state that our 
initial conditions do not refere to present cosmic time (usually labeled as t = to), and so we will not use the subindex 

a 5? 
■ 

The spherically symmetric LTB metrics (5) are contained within a larger class of more general metrics (the Szekeres- 
Szafron metrics j30|, admitting in general no isometries. The integration of the field equations for (5), given a 

source (3) satisfying (2), is examined in the next section. For the case of more general Szekeres-Szafron metrics, see 
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IV. INTEGRATION OF THE FIELD EQUATIONS. 



Einstein's field equations for (5) and (3) are 



up = — 



Y (f 2 + Fc 2 } 



Y 2 Y' 



-G* j 



(9a) 



up 



Y Y 2 + Fc 



2Y 2 Y 



3Y 2 Y' 



1(2G% + GV) , 



(96) 



kP 



Y 



Y (t 2 • Fc 2 ) -+ 2 V " : V 



(9c) 



where k = 8ttG/c 2 . Imposing on (9a) and (9b) the equation of state (2), using (8) and integrating with respect to r 
yields the following constraint 



2Y(Y 2 + Fc 2 ) + Y 2 Y -Kmc 2 I n^Y 2 Y^dr = A(t) , 



7' 



(10) 



where X(t) is an arbitrary integration function. It is important to remark that (10) follows only from (9a) and (9b) 
without involving (9c), ie it was not necessary to make any assumption regarding the form of P in order to obtain 
(10). A second integration of the field equations necessarily requires setting X(t) = in (10), leading to 



M + W { y 



Fc 2 , 



(11) 



where 



M 



( TO ) 2 (m) 

Pi =mc n, 



(12a) 



W 



= J p[ r) Y 2 Y>dr, pP^SnT^Ti 



(126) 



so that p\ m \ p[ r ^ respectively define the initial densities of the non-relativistic and rclativistic components of the 
mixture. 

In the remaining of the paper we restrict ourselves to F = 0, similar to the choice of spacelikc sections of zero 
curvature in FLRW geometry, leaving the case F ^ for a future analysis. An explicit integral of (11) in this case is 
given by 



-^{t - U) = V^+~e{y - 2e) - VT+~e{l - 2e) , 



where 



p = 



kM 
Y?' 



e = 



W 

M' 



y = 



Y 

y: 



(13a) 



(136) 



It is possible to invert (13a), thus obtaining y = y(t,r) as a complicated, but closed analytic form, where the r 
dependence is contained in the functions /i, e appearing in (13b). However it turns out to be more convenient to 
use (11) and (13) to simplify the field equations and radial gradients of y in order to express all state and geometric 
variables in terms of y and suitable initial value functions related to those of (12) and (13b). 



•5 



V. THE STATE VARIABLES. 



From (8) and (11)-(13) it is possible to obtain the state variables 

n (m) ) n (r) ) T,p {m \ p {r \ p, P. However, before 

doing so it is useful to define the averaged initial densities 



Pi I — y? Y 3 ' 



Ar) 



Jp^d(Y 3 )_3W 



y. 3 ' 



(14) 



averaged over the volume Y 3 . Since the solutions allow for an arbitrary re-scaling of the radial coordinate, without 
loss of generality we can select Yj = rR i7 where Hi is a characteristic constant length scale. Therefore, the volume 
Y 3 , evaluated from the symmetry center, r = 0, to an arbitrary fluid layer r, can be characterized invariantly as the 
volume of the orbits of the rotation group SO (3) in the hypersurface t = t,. In the Newtonian limit, the distance Yi 
becomes the radius of the circular keplerian orbit in the field of (5). 

Together with the averaged initial densities, we shall define the quantities A- m \ A- r ^ given by 



(m) 



(m)i 



T 3 dr 



Pi 



(m) 



3M 



A™) 



- 1 



Am) 



1 + A. 



(m) 



(15a) 



A 



(r) 



JbPYYfdr _ p 



(r) 



3W 



Ar) 



(r) I (r) 
Pi = {Pi 



1 + A 



(156) 



whose interpretation as effective initial density contrasts is discussed in the following section. Using (14) and (15) we 
can re-write p, e in (13b) as 



3Y°* I 3(1 + Ar») : 



Pi 



(r) 



Pt" 1 + 4 



Am) 



1 + A 



(r) 



(16) 



The state variables n^ m \ n^ r \ p, p, P now follow by inserting (11) into (8) and (9), while T is obtained with the help 
of (2). This yields the following forms 



Am) 



(m) 



y 3 r' 



Ar) 



(r) 
n i 

y 3 r' 



(17a) 



T 

T= — 

y 



(17b) 



p = p^ + p& 



yi yi 



1 



(17c) 



P 3y 4 r ' 



(17d) 



P = 



(r) ^ 
Pi g 

6y4 r ' 



(17e) 



where the functions T, and $ are given by 
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Y'/Y 



# = 1 



(i-r) 



$ = 1 + 



(1 - 4T) 



(18) 



W * * ' 3(1 + A^)' " ' ' 3(1 + A< r) ) 

The solutions characterized by (2), (4), (5)-(18) become determinate once T above is obtained in terms of y from 
L3) for g 
(18) into 



(13) for given initial value functions p\ m \ p^ (re-expressed in terms of the quantities A-"*' , A- r) , e). This transforms 



r = 1 + 3AA-" 1 ' + 3BAj P) , 



* = 1 



AA^+BAV 



1 + A 



(r) ! 



(19a) 
(196) 



where 



$ = 



-4AA\ m) + (1 - 45) Af 



1 + A 



(r) 



(19c) 



1 

V 



,r I-// 8e 2 



(1 - 4e - 8e 2 ) 



(20a) 



2/ 



y + 2e- (l + 2e) 



(206) 



with A< m) , A| r) given by (15). The kinematic parameters a, 9 follow by inserting (11) with y = Y/Yi and (19a) into 
(6a) and (6b) 



a = — - 





Ay A T 


} +B,yA, 


r)" 


y 2 


1 + 3 


(AA[ m) 


+ BA<r) 





(21) 



e VmV?/ + e 


1 + (3A + Af ' + (3B + yB. y ) Af ' 


3 y 2 


l + S^A^+BA^)" 





(22) 



where i? y are the derivatives of A, B in (20) with respect to y. Given a set of initial conditions specified by 

e, A/" , A[ , equations (17) and (19)-(22) provide fuly determined forms of the state and geometric variables as 
functions of y and the chosen initial conditions. 

The solutions presented so far contain a FLRW particular case, obtained by setting in (11) and (12) n-" 1 ' = n[ m \ 

aif^ = nf^ and Tj = fj, where n\ m \ nf\ Ti are arbitrary positive constants. Under this parameter specialization, 
(13) holds with Y = R(t)f(r) (so that y = R(t)) and (5) becomes a FLRW metric. This leads to 




r = * = i, $ = o, 



where p\ = mc 2 nl , pf — 3n\ r 'k B Ti, so that T = T(t), p = p(t), p = p(t) and P = 0, with (3) becoming a perfect 
fluid tensor where p and p satisfy (2). The FLRW limit can also be characterized by A| ' = A\ ' = 0, and so, from 
(11), equations (21) and (22) become: a = and 6/3 = y/y — R/R. Another limit is that of LTB dust solutions, 
obtained by setting T l = in (12b) and (11), so that (17) becomes T = p = P = and p = mc 2 n' m) . 
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VI. DENSITY CONTRASTS AND REGULARITY CONDITIONS. 



Since the radial dependence of all state and geometric variables is sensitive to ' and ' defined in (15), it 
is important to provide an interpretation for these quantities. From (14) and (15), it is evident that A^ ' and A^ 

are effective "gauges" of the deviation of p± , p'p from their volume averages for every closed interval in the range 
of the integration variable r along the initial hypersurface t — ti. The signs of these quantities characterize initial 
density profiles, with "density lumps" as these densities decrease {[p^]' < 0, [p^]' < 0) or "density voids" as they 

increase ([p\ m) ]' > 0, [p\ r) }' > 0). Also, with the help of Rolle's theorem applied to (14) we find that A i and A i , 
as adimensional functions of r, are constrained by the maximal density contrasts in terms of 

( m )max 
| A- | < P \ . , - 1 , 

r % 

Wmax 



A < 



„( r )min 

ri 



(m) (r) 

where the superindices "max" and "min" respectively indicate the maximal and minimal values of p\ , p\ in any 
interval < r along the hypersurface t = ti. Small initial density contrasts obviously imply 

p (m)max ^ p Mmin / ^ |Af , |«l, (23a) 



Mmax^ (r)min ) (r) M /(r)\ ^ | A f|«i, (236) 



3 fX m > r.( m ) mc z mc z 

Pi 11 i 

allowing us to consider a formal analogy between A^ ' and A^ and energy density "exact" initial perturbations. 



This is further reinforced from the definitions in (15), and by remarking that the FLRW "background" follows by 

"turning the perturbations off", that is, setting: A^ = A[ = 0. 
An important restriction that the solutions must satisfy is the following regularity condition 

Y'/Y 



Y'/Yi 



> 0, (24) 



which prevents negative densities n {r \ n(™\ as well as the occurrence of a shell crossing singularity p5j. This 
singularity is characterized by unphysical behavior because Y appears in the denominator of equations (17a), (17c), 
(17d) and (17e), but does not appear in (17b). Therefore, if T = 0, the densities, pressure and viscous pressure 
diverge with T finite (in general), a totally unacceptable situation that can be avoided by considering only the range 
of evolution of the models to spacetime seccions with t > U satisfying (24). The fulfilment of (24) depends on the 
functions A, B and on the magnitudes of the initial density contrasts gauged by and A^ m) . This will be examined 
further ahead together with the conditions for thermodynamical consistency. 

VII. THERMODYNAMICAL CONSISTENCY 

The models derived and presented in the previous sections must be compatible with a suitable thermodynamical 
formalism. For this purpose, it is advisable to leave aside the "conventional" theory of irreversible thermodynamics 
p6|, fl37| , whose transport equations are unphysical as they violate relativistic causality of the dissipative signals as 
well as stability of the equilibrium states (seee.g. p8), |M and p3]). We shall consider instead "Extended Irreversible 
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Thermodynamics" (EIT) poj , ]2l| ], p2| , [ p3| , a theory free of such serious drawbacks and so a more adequate 
theoretical framework for the models. According to EIT, the entropy of a system away from thermodynamical 
equilibrium depends not only on the "conserved" (equilibrium) variables (i.e. particle number densities, energy 
density and so on), but also into the non-equilibrium fluxes (i.e. heat flux, and bulk and shear dissipative stresses). 
This theory is supported by Kinetic Theory of gases, Information Theory and by the Theory of Hydrodynamical 
Fluctuation - see |g3[ for a detailed description. When shear viscosity is the only dissipative agent, the corresponding 
generalized entropy s of radiation plus matter obeying the usual balance law with non-negative divergence, up to 
second order in II ab takes the form 



8 = 



nT 



run 



ab 



(snu a ) ;a > 0, 



(25) 



where is obtained from the integration of the equilibrium Gibbs equation, n = n^ r ' +?i/ m ) and a is a phenomeno- 
logical coefficient to be specified later. The evolution of the viscous pressure is, in turn, governed by the transport 
equation 



rtl cd h c a h d b + n 



ab 



X rr ( T < 

— 1 11 — — u 
2 '\Tr, 



2r]a ab = 0. 



(26) 



where i] and r are the coefficient of shear viscosity and the relaxation time of shear viscosity, respectively. The former 
as well as other related quantities can be obtained by a variety of means Q including Kinetic Theory, Statistical 
Mechanics or both |^5| , 26 1 , |27]] , |^8| . The relaxation time, r , is related to and larger than the mean collision time 
between particles and it may, in principle, be estimated by collision integrals provided the interaction potential is 
known. As a physical reference to infere the form these coefficients might take, consider the "radiative gas" , with p, p 
satisfying (1) (or the approximation (2)). For the radiative gas the forms of i], a in terms of the relaxation time of 
the dissipative process r are 



4 (r) 



2r?, 



8p« 



(27) 



where p^ r > is either n^k B T or a^T 4 and the subscript "(rg)" emphasizes that these quantities are specific to the 
radiative gas. 

We verify now the compatibility of (13)-(22) with (25)-(27). Integrating the equilibrium Gibbs equation and 
substituting (27) into (25), we obtain 



1 fc D ln 



3n 



(r) 



nr (T 



15k ( P 
P 



4a n I 3 



3n 



r) 



k a In 





15fc 




r 


~ "32~ 





(28) 



where the approximation: n = n ' +n w n r was used and the initial value of has been set to be the equilibrium 
entropy per photon. Equation (28) reflects the fact that we are neglecting the contribution of the entropy due to 
non-relativistic particles, a justified approximation since the latter are much less abundant than the photons. 

Ideally, the transport equation (26) should be satisfied for T) having the form (27), associated to the radiative gas, 
and the relaxation time, r, given by collision integrals obtained from Kinetic Theory. However, as mentioned in the 
introduction, and in order to obtain exact expressions for all thermodynamical parameters, we will assume r\ given 
by (27) and deduce r from the fulfilment of (26). This yields 



1 + A 



3 + 4Aj r ) + gr 



~ 256 



(29) 
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While there is no need to justify r\ given by (27), this form of r is acceptable as long as (29) satisfies the requirement of 
a relaxation parameter: it must be a positive quantity and must comply with a positive entropy production law s > 0. 
It is desirable that r should somehow relate or approach its definition as a collision integral and that its behavior be 
qualitatively analogous to a suitable mean collision time, therefore it should be an increasing (decreasing) function if 
the fluid is expanding (collapsing). Evaluating s from (28) and comparing with (29), we obtain the following relation 
between s and r 



15k a fP\ 2 15fc„ /$ x 2 



4r \p J 16t V* 



(30) 



consistent with the general relation |^] (nsu a )- a — H a bH ab / (2rjnT) associated with (25) and (26). As a consequence 
of (30), s > and r > imply each other. Also, from the form of (29), necessary and sufficient conditions for positive 
t, s, p, T can be given by 

* > 0, (31a) 

cr$ < 0, (316) 

while the condition ensuring concavity and stability of s can be phrased for an expanding fluid configuration as the 
requirement that s decreases for increasing r (f > s < 0). From (30)-(31), this follows as 



f > 0, 



s 2aT {Pi 



s 3*$ >) 

ri 



0) 
Pi 



Q O) 



- < 0. (31c) 

T 



So that, if (24), (31a) and (31b) hold, then (31c) reduces to f > 0. 

Since r is a thermodynamic relaxation parameter, it is important to compare it with another natural timescale of 
the models: the Hubble expansion time defined by 

t H = |. (32a) 

where follows from (22). Such comparison should provide an insight into the timescales associated with the mixture 
interaction and decoupling. However, strictly speaking, the criterion for interaction and decoupling in cosmological 
gas mixtures is not given by comparing r and t H but by comparing the latter with the timescales associated with the 
various reaction rates of the radiative processes involved, particularly the photon mean collision time t obtained from 
Thomson scattering ||, Hence, we can consider this relaxation time as approximately gauging the interactivity 
of the matter mixture by demanding that for a range of the evolution of the mixture, approximating its interactive 
range, we must have 

t <t H (326) 
while, as the mixture evolves and the components decouple, eventually 

r>t H . (32c) 

Since r in (27) and (29) must behave qualitatively similar to t [^9|, the comparison in (32b-c) should yield qualitatively 
analogous results as a similar comparison between t H and t . The temperature associated with the passage from (32b) 
to (32c) (obtained from (17b)) should approximate the decoupling temperature obtained by the condition t = t H . 
This point is examined in section XI. 

Equations (31)-(32), together with the regularity condition (24), provide the necessary and sufficient conditions 
for a theoretically consistent thermodynamical description of the solutions within the framework of EIT and Kinetic 
Theory applied to the radiative gas. We examine the effect of these conditions in the following section. 
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VIII. THERMODYNAMIC ALLY CONSISTENT MODELS 



A. Conditions (24) and (31a) 

From (19a) and (19b), the fulfilment of (24) and (31a) is equivalent to the following condition 

-\<AA\ m) +B/^ r) <1 + A l r\ 



(33) 



From (20), since the functions A and B diverge as y — > there are necessarily values of y for which (33) is violated. 
However, from (19a-b) we have: = ^ = 1, and the range of y we are interested is y > 1. As shown by figures la 
and lb, displaying the implicit plots of T = and ^ = for an initial temperature Ti ss 10 6 K (from (23c), e sa 10 3 ), 

there are no zeros of T and * for y > 1, |A- | < 1 and |A- | < 1. Therefore, (33) holds in the range of interest for a 
large class of initial conditions. 

B. Conditions (31b), (31c) and (32) 

For the examination of these conditions we will assume that (33) holds for y > 1, | A i " 1 | < 1 and \ < 1 (see 
figures la and lb). Then, with the help of (19)-(22), condition (31b) is equivalent to 



C = AAA . 



(m) 



+ - l)B,y 



(m) (r) 



+ [4AB >v + (4B-l)A tV ]A i A 4 < 0. 



An insight into this expression follows by looking at its initial value 



Ci = 



-A 



2(1 + e) 



Af+eAr 



and its asymptotical behavior as y 3> 1 



(34a) 



(346) 



Ac 



C 



A 



V 



+ 



4AA 



+ 



35 eAA 



(») 



VTT^y 5 / 2 9 VTTiy 7 / 2 



Ai 



A 



00 



+ 24A 2 A^A 



24(1 + e) 



where 



and 



+ 



64 e 2 AA 



Ai 



A. 



24A 2 A, (r) A( s ) +9 



A 



63 yi+ly 9 / 2 18 



(l + e)2/ 4 



55 e 3 AA, (s) 



8 ^/TTiy 11 / 2 
A (s) 



|A (r) -A <m) , 

4 » "» > 



(34c) 
(35) 



^ = A - m> A 2 — 3e(l + 2e) A , 

Ai = l-8e + 2e 3 + 2e 4 , A 2 = 8e 2 +4e-l. 

As shown by (34c), the quantity defined in (35) plays a fundamental role with regards to the fulfilment of (31b) 
and (31c). This suggest classifying initial conditions in terms of A^ s \ We shall examine (31b), (31c) and (32) for 



each case A^ = and A- SJ ^= separately. 
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1. The case A\ s) = 

^From (34c), the condition A- s ^ = is necessary and sufficient for having C < asymptotically. In fact, substituting 
A| s ' = into (34a), leads to the following dramatic simplification of C 



C 



A Ml2 



(3y + 4e) 

(36) 



8(1 + e) y 5 



so that A^ — is a sufficient condition for the fulfilment of (24), (31a) and (31b) along the range y > 1, for |A| r> | < 1 
and lAfl < 1. 

Condition (31c) is also satisfied, since r increases monotonously along the fluid worldlines, behaving asymptotically 

as: r w y 3 ^ 2 - Hence, those models whose initial conditions satisfy A- s -* = can be characterized as the subclass of 
models complying with the conditions of thermodynamical consistency in the asymptotic range of y. This point is 
consistent with the fact that A^ = is a necessary and sufficient condition for having \P\/p — |$|/(2\I>) — > and 
s — > as y — > oo, so that the fluid layers evolve towards an asymptotic equilibrium state. However, using t H =3/9 
calculated by inserting A^ = into (22), we have r/t H < 1 for all the evolution of the fluid, thus failing to comply 
with (32b-c). This is illustrated by figure 2, and implies that the relaxation time r cannot be associated with a 
radiation matter mixture whose components interact and then decouple. 



2. The case A (s) / 

If A( s ) =/= 0, condition (31b) cannot hold along the full range of y (because of (34b)), but might hold along a 
restricted range of physical interest 1 < y < y» for which the mixture could be in the interactive stage. This situation 
is not incompatible with the thermodynamical arguments of the previous section, since the phenomenology of the 
radiative gas model strictly applies if the mixture components interact. The fact that condition (31b) can hold for 
1 < V < V* follows from evaluating the sign of C, given by (34a), a quantity that can be negative (so that Tj > 0) for 
a wide range of acceptable situations (for example, if A^ ' and a\ ' have the same sign). However, as y increases, r 
either changes sign or diverges positively, depending on the zeroes of $ and a. Since we are assuming that (24) and 
(31a) hold, the sign of r, as given by (29), depends only on the quotient $/ct, and so the behavior of r is strongly 
related to the signs and zeroes of these functions. If > and, as y increases, there is a zero of $ for a ^ 0, then 
r passes from positive to negative, but if the zero of a appears first, then r diverges positively. A zero of <I> (with 
cr / 0) might be compatible with (31b) (r > 0), but violates (31c) and so is unacceptable. However, a zero of a (with 
$ 7^ 0) is acceptable, since r would diverge positively (as y — > y*) and so would be a positive and increasing function 
along the range 1 < y < y*. In order to verify if this type of evolution is possible, it is necessary to gather information 
on the zeros of $ and a. Since these expressions are cumbersome, it is convenient to examine their zeros graphically, 
and so we have plotted in figures 3a and 3b the solutions of the implicit equations $ = and a = 0, in terms of A^ , 
A i and log 10 (y), while figure 4 displays the sectors in the plane A i , A i where $ = and o = occur. As these 
figures reveal, sufficient conditions for the desired evolution are given by 

A| s) > for Af > (37a) 

A| s) < for Af < 0. (376) 

Therefore, if A| s ' ^ 0, conditions (37) are sufficient for the fulfillment of (31b) and (31c) along the range 1 < y < y», 
where j/» is a zero of a. Conditions (32) are satisfied under the restrictions (37). This will be discussed further ahead 
in section X. 



IX. INITIAL CONDITIONS AND EXACT INITIAL PERTURBATIONS. 

The behavior of the quantity A^ given by (35) determines the set of initial conditions that characterize the 
thermodynamical consistency of the models. Since this quantity plays such an important role, it must be related to 
a physically significant property along the initial hypersurface t — ti. From (35), using (12), (14) and (15), we obtain 
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A« = ?A, W 



A 



1 [log(^ 3 / 4 /M)] ' _ dflogE, 



pog(K 3 )]' 



d [log if]' 



(38a) 



where 



3/4 



(m) 



(386) 



In order to provide an interpretation for (38), we remark (from (17), (19) and (28)) that the entropy per photon along 
the initial hypersurface 



4a D T? 15 



3?) 



(r) 



32 



2 Aa^Tf 15, 
— k c 



3n 



(r) 



32 



A 



is very close to its equilibrium value, since its off equilibrium correction is proportional to [A 4 ] 2 , a very small quantity 



since A^ is already assumed to be small. Therefore, using the approximations associated with small density contrasts 
(equations (23a-c)), we obtain 



( r ) ~ O T^4 



E oc 



A™) 



,(m) 



Am) 



(39) 



implying that Hi is proportional to the ratio of the averages of photon entropy and baryon number density, which 

(for small density contrasts) is roughly equivalent to the averaged ratio of these quantities. Hence, condition = 
roughly means an initial hypersurface with constant averages of photon entropy per baryon, while condition (37) 
roughly means (see equations (20)) that the sign of the spatial gradient of the average of photon entropy per baryon 
must agree with the sign of the gradient of the initial photon energy density: as r increases along t = ij, it must 
increase for a density void Qfi^ ]' > 0) and decrease for a density lump ([p-^]' < 0). 

The assumption of small density contrasts leads to a natural comparison with the Theory of Perturbations of a 
FLRW background, in the isochronous gauge and considering a mixture of radiation and non-relativistic matter, see 
f§! (§L 07 |§j @; @- Under these assumptions, matter and radiation densities are given by 



p (rn) _ -(m) fj + S (rn) 



(40a) 



p (r) = - p (r) 



1 + 5^ 



(406) 



where p {m \ are the respective densities in a FLRW background and < 1, < 1 are the perturbations. 
The gauge invariant quantity 



(41a) 



which formally resembles (35), defines the fluctuations of photon entropy per baryon and leads to the classification of 
perturbations as 



"Adiabatic" 



5^ s) = 0, 



"Isocurvature" 



(m) 



3p 

4p(r) 



(416) 
(41c) 
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In order to establish a comparison with Perturbation Theory, consider defining the initial densities along the lines of 
(40) 



J«0 



Am) 



1 + 6, 



(m) 



Pi 



(m) 



r(™) 



l + (6. 



r(m) 



(42a) 



(r) _(r) 
Pi = Pi 



1 + 5 



(■'•) 



(r)\ _(r) 
Pi > = Pi 



;(r) 



(426) 



where p) , p) are now the constant values of the initial densities in a FLRW background and | <^ m) |< 1, | 5^ |< 1 
are exact initial perturbations. Inserting (42) into (38) leads to 



Am) 



; (m) 



tf 5 



4 i 1 



dK 3 



(42c) 



where Si oc [p| r ''] 3 ^ 4 and n| r ' are the photon entropy and baryon number density of the FLRW background in the 
hypersurface t = U. The rhs of (42c) illustrates that, under the assumptions (42a-b) of a perturbative treatment, 
approximately reduces to the radial gradient of 5^ evaluated in t = ij. However, the comparison with perturbations 
must be handled with caution, since in the case of (35) and (36)-(39) we are dealing with quantities and relations 
strictly defined in the initial hypersurface, and so determining initial conditions in terms of averages of initial value 
functions. The Theory of Perturbations, on the other hand, traces the evolution of 5^ m \ for all t. In particular, 
we must be specially careful if we borrow the terminology of perturbations in (41) to characterize the cases A^ = 

and A^ s) ^ 0, since in our case a relation like A^ = specifies initial conditions and will not be satisfied (in general) 
for t > ti. Another delicate point refers to the current use of the term "adiabatic" for the case (41b), meaning 
perturbations that conserve photon entropy. Strictly speaking, these perturbations should be denoted as "isentropic" 
or "reversible", since an adiabatic process need not be isentropic (or reversible). It is quite possible that this conceptual 
vagueness follows from the fact that most papers in Perturbation Theory, either assume thermal equilibrium, or have 
incorporated dissipative processes without the necessary rigour. Hopefully, recent work [ |l7| on these lines might be 
helpful to clarify these issues. Raising this point is relevant because the models presented in this paper assume a fluid 
evolving along adiabatic but irreversible (non-isentropic) processes, therefore initial conditions like A^ = (formally 
analogous to (41b)) do not imply a conserved photon entropy, but a roughly constant average of photon entropy at 
t = ti, but not (in general) for t > U. Never the less, for the sake of mantaining continuity with currently established 

terminology, we shall use the formal analogy between (41) and (42) in the examination of the cases A^ s) = and 

(s) 

Aj 0. Developing further the comparison with perturbations on a FLRW background is relevant and interesting. 
However, such a task requires a comprehensive and detailed elaboration and so will be carried on elsewhere. 



X. THE CASE A { j S) ± 0, AN INTERACTIVE MIXTURE. 

From section VIII we know that if (37a-b) hold, then r diverges positively as y — > where is a zero of a. We 
still need to know if rjt n < 1 along the range 1 < y < y, where the set of values y = y < y* are characterized 
by r/t H = 1, approximately marking the decoupling of matter and radiation. A sufficient condition for this type of 
evolution follows from evaluating [r/t H ]i at the initial hypersurface 



5A« 


1 + A^ 


"(3 + 4e)A l (r) -4A J (s) +8(l 


+ s)_ 


_ 10(1 + e) _ 5 


"(3 + 4e)A l (r) -4A l (s) (3 + 4)" 




36 + 47A l (r) + 16 




2~ 


" 9(3 + 4e) ~ 18 



<1. 



where we have eliminated A^ m) from (35) and assumed A^ <§; 1, A^ r) <§; 1 and e « lO -3 ^ « 10 3 . Since r diverges 
as y — > but O remains finite in this limit, the ratio T/t H , initially smaller, necessarily becomes larger than unity 
for 1 < y < y < y*. 
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The currently accepted value of matter and radiation decoupling is T D w 4 x 10 3 K. Assuming w 10 6 K, the set 
of values y = y D associated with this temperature follow from (17b) by solving for y = y D the equation 

T D «4xl0 3 = ^ (2/£) ,A( s) ,Ar), (43) 

where *f? is given by (19b). Assuming small density contrasts, figure 5 illustrates that this value of T D is closely 
associated with y D w 10 2 4 . Also, it is evident from figure 3b, that having a zero of a for values 10 2 < y» < 10 3 

requires a very small deviation from = 0. This is illustrated by the approximated sketch of the level curve 
y = y D = 10 2 - 4 that appears in figure 4, implying that y, and y D lie in a a very narrow sector of the plane 
A^ , very close to the line A^ = 0. Therefore, for these values we must have y* ~ y ~ y D - However, it is not 

(s) 

clear from figures 3b and 4 how small the deviation from A^ = should be. Hence, we have ploted in figure 6 the 
implicit equation a — 0, showing a more precise relation between the orders of magnitude that the occurrence of the 
zero of a for y D w y» ~ 10 24 implies for |A^ r) |, |A^|. 

Another constraint the models have to comply with is the observat ional bounds , JgJ on the anisotropy of the cosmic 
microwave radiation (CMR), given by the maximal photon temperature contrast [<5T/T] D sa 10~ 5 . For y w 10 2 ' 4 and 
Ti w 10 6 , so that e sa 10 3 , we have: A w 1.1 x 10 _1 and B sa 2.2 x 10 _1 , were A and B are given by (20). Also, from 
(17b) and (19b), the deviation from the equilibrium FLRW form T eq = Ti/y at y = y D is given by 

._ 1.1 x 10"^"° + 2.2 x lO^Af fAA . 
*-l « l - — i— , (44) 

indicating that compliance with the maximal observed temperature contrast of the CMR constrains the maximal values 
of A^™', A,| r) to about 10~ 4 . Therefore, from figure 6, the corresponding variation range of |A^ S '| is |A^| < 10~ 8 , 
and so compatibility with acceptable values of T D and the CMR anisotropy implies |A^ S ' | sa \Al ' | 2 <C |A^ ' |. 

Under the analogy between A^ ' , A\ ' with perturbations in a FLRW background, the case A^' ^ seems to 
correspond to isocurvature initial perturbations. However, since p^ m > <C before decoupling, the latter are char- 
acterized (see equation (41c)) by 5^ « and 5^ <C 5^ s \ and so, following the analogy between (41) and (42) 

(s) 

would disqualify A \ ^ as comparable to an isocurvature initial perturbation. In fact, the maximal bounds on 

A^ ' and A^ obtained in the previous parragraph yield exactly the opposite behavior. Borrowing the terminology 
of Perturbation Theory, the analogy between (41) and (42) would characterize "isocurvature" initial conditions by 

A/ =0. From figures 3 and 4, this condition is not incompatible with (31) and (32), but yields a decoupling surface 
y D ~ 10 with a decoupling temperature much larger than the accepted value. Hence, the acceptable values of the 

( s) 

case AJ / cannot be associated with this type of initial perturbations, but since perturbations are in general 

combinations of adiabatic and isocurvature components and we have A[ s ' <C A i ' and A^ s ' <C A i , a more accurate 

analogy for A^ ' , A^ ' is that of "quasi adiabatic" initial perturbations. Following this analogy, and bearing in mind 

that A| m) and Al ' are defined for an initial hypersurface characterized by TJ ~ 10 6 K. (still outside the horizon for 

perturbations of all wave numbers), the bounds on the magnitude of A,-*" 1 can be related to bounds on the amplitude 
of "nearly adiabatic" fluctuations of photon entropy per baryon generated in the earlier, more primordial, inflationary 
era 

B, 0, @, 01, II- Primordial entropy fluctuations have been examined in connection with the creation of 
axions in an inflationary scenario j^], JlO| |. However, these are isocurvature perturbations, and so might not be related 

to entropy fluctuations associated with A^ ^ 0. The study of entropy fluctuations in inflationary models is still a 
highly speculative topic fl34| , and its connection with initial conditions in the models under consideration should be 
an exciting subject to examine in a future paper. 

XI. SAHA'S EQUATION AND THE JEANS MASS 

As mentioned earlier, it is important to compare r and t H with the timescale characteristic of the interaction rate of 
the photons and electrons for Compton and Thomson scattering (the dominant radiative processes in the temperature 
range of interest) 
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K = ' *c = -rVi T , (45) 

carn e k B l 

where cry ~ 6.65 x 10~ 25 cm 2 is the Thomson scattering cross section, m e is the electron mass and n e is the number 
density of free electrons, a quantity obtained from Saha's equation 
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(46) 
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where X e is the fractional ionization, h is Planck's constant, ng is the number density of baryons and B a 
is the binding energy of the hydrogen atom. Combining (45) and (46) we obtain 



13.6 ev 



1 



2ca n 



4h 3 



(m) \ 1/2' 

n exp(B /k B T) \ 



(27TTO e k B T) 



3/2 



(47) 



The recombination process is characterized by X e ss 1/10 in (46), so that most free electrons have combined with 
protons into neutral atoms, while the decoupling of matter and radiation strictly follows from the condition: t = t H , 
a condition analogous to t — t H and leading to the "decoupling temperature" from (45) and (47) with the help of 
(17b) and (22). Since t must be smaller but qualitatively similar to r (increasing as r increases j29|), a comparison 
between r and t H for near decoupling temperatures should be qualitatively analogous to that between t and t H . 
The timescales t and t c can be easily obtained from (45), (47) and (22) and compared with t H , t for a set of 
parameters characteristic of thermodynamically consistent models discussed in the previous section. As shown by 



figure 7, displaying the ratios of r, i 7 , t c to t H for A 



(r) 



10" 4 and 10" 8 < A, (s) < 1CT 4 . For A 



10 and for near 



decoupling temperatures (10 < T < 10 K. or equivalently 10 < y < 10 ) we have t smaller than r (but of the 
order of magnitude) and overtaking t H along a set of values of y that closely match the accepted value T D w 4 x 10 3 K. 
This is consistent with the estimation: A^ w 10 -8 obtained in the previous section (see figure 6). On the other hand, 
for higher temperatures (T sa 10 6 K. or equivalently y ~ 1), the relaxation parameter r is of the order of magnitud 
but greater than the Compton timescale, the dominant radiative process at these temperatures. It is also possible to 
show that X e sw 1/10 in (46) leads to the accepted values of the recombination temperature. 



To end the discussion, we compute the Jeans mass associated to the initial conditions of the case A 
mass is given by §, §, 0, §, | 



^ 0. This 
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(48) 



where p, p, are given by (2), (7), (13) and (14), Xi = PT/Pi' * an d T follow from (19) and C s is the speed of 
sound, which for the equation of state (2), has the form 



C =^r 



3p 



(m) 



pi™) = mc 2 ni m 



\ p( r > = 3n< r >fc B T. 



Evaluating (48) for y = y D w 10 2 4 , e w l/xi ~ 10 3 and p { [ ] w a B T" 4 w 7.5 x 10 9 ergs/cm 3 , yields M, w lU'-gm, or 
approximately 10 16 solar masses. This value coincides with the Jeans mass obtained for baryon dominated perturbative 
models as decoupling is approached in the radiative era. 
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XII. CONCLUSIONS 



We have derived a new class of exact solutions of Einstein's equations providing a physically plausible hydrody- 
namic description of a radiation-matter (photon-baryon) interacting mixture, evolving along adiabatic but irreversible 
thermodynamic processes. The conditions for these models to be consistent with the transport equation and entropy 
balance law of EIT (when shear viscosity is the only dissipative agent) have been provided explicitly, and their effect 
on initial conditions have been given in detail, briefly discussing the analogy of these conditions with gauge invariant 
initial perturbations in the isochronous gauge. As far as we are aware, and in spite of their limitations mentioned 
in the introduction, we believe these models are the first example in the literature of a self consistent hydrodynam- 
ical approach to matter-radiation mixtures that: (a) is based on inhomogeneous exact solutions of Einstein's field 
equations, and (b) is thermodynamically consistent. We believe these models can be a useful theoretical tool in the 
study of cosmological matter sources, providing a needed alternative and complement to the usual approaches based 
on perturbations or on numerical methods. 

The solutions have an enormous potential as models in applications of astrophysical and cosmological interest. 
Consider, for example, the following possibilities: (a) Structure formation in the acoustic phase. There is a large 
body of literature on the study of acoustic perturbations in relation to the Jeans mass of surviving cosmological 
condensations. Equations of state analogous to (2) are often suggested in this context §, @, @, @, §, @. 
Since practically all work on this topic has been carried on with perturbations on a FLRW background, the exact 
solutions derived and presented here may be viewed as an alternative treatment for this problem, (b) Comparison with 
Perturbation Theory. The models presented in this paper are based on exact solutions of Einstein's field equations, 
but their initial conditions and evolution can be adapted for a description of "exact spherical perturbations" on a 
FLRW background, ft would be extremely interesting, not only to compare the results of this approach with those of 
a perturbative treatment, but to provide a physically plausible theoretical framework to examine carefully how much 
information is lost in the non- linear regime that falls beyond the scope of linear perturbations. We have studied in this 
paper only the case F = in (11), thus restricting the evolution to the "growing mode" since all fluid layers expand 
monotonously. The study of the more general case, where F(r) in (11) is an arbitrary function that could change sign, 
would allow a comparison with perturbations that include also a "decaying mode" related to condensation and collapse 
of cosmological inhomogeneities. (c) Inhomogeneity and irreversibility in primordial density perturbations. The initial 
conditions of the models with ^ are set for a hypersurface with temperature Tj w 10 6 . These initial conditions 
can be considered the end product of processes characteristic of previous cosmological history, and so the estimated 
value |A^| ~ 10~ 8 , related to the spatial variation of photon entropy fluctuations, can be used as a constraint on 
the effects of inhomogeneity on primordial entropy fluctuations that might be predicted by inflationary models at 
earlier cosmic time. Also, the deviation from equilibrium in the initial hypersurface (proportional to [A^ '] 2 « |A] S '|) 
might be helpful to understand the irreversibility associated with the physical processes involved in the generation of 
primordial perturbations These and other applications are worth to be undertaken in future research efforts. 
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FIG. 1. The equations T = and = 0. Implicit plots of the solution of V = (figure la) and = (figure lb), in 
terms of A; ' , A; ' and log 10 (jy), with V, given by (19). The grid marks the initial surface y = 1, and is not intersected by 
the surfaces r = 0, * = 0, illustrating that V and ^ have no zeroes in the evolution range y > 1 for density contrasts bound 
by — 1 < a[ ' < 1 and — 1 < A^ ' < 1. We have considered /n^ m ' « 10 9 , m to be the mass of a proton and T; ~ 10 6 K, 
hence we have made the approximation (23c): e ?a 10~ 3 Ti w 10 3 . These (and the remaining) plots were obtained with the help 
of the symbolic computing package Maple V[39]. 



FIG. 2. The ratio of relaxation vs Hubble times for the case A^ s) = 0. This plot displays the ratio r/t H , as a function 
of log 10 (y) and log 10 (|A; ' |) (denoted as "log(Dr)"), where r is the relaxation time and t H — 3/0 is the Hubble expansion 
time, given by (22) and (29) under the condition A^ = 0. This ratio is less than unity for all y, and |A ; ' | < 1 thus violating 
(32b) and (32c). 

FIG. 3. The equations <3> = and a = 0. Implicit plots of the solution of $ = (figure 3a) and a = (figure 3b), in terms 
of a[ , A, and log 10 (y). The "wall" where the surfaces a — and $ = occur for values greater than y — 10 2 corresponds 
to the line jA[ ' — a[ ' (or A* s ' = 0). The plots also illustrate that the values of A^ ' , A^ ' where a zero of a occurs are 
clearly distinct from those associated with a zero of leading to conditions (37). The vertical height y — 10 2 is displayed in 
figures 3a and 3b, illustrating that the values of Aj ' , a\ ' with 10 2 < y* < 10 3 are very close to A^ = 0. 



18 



FIG. 4. Allowed values of A- , A; for a = 0. The fi gure displays the plane A, , A, associated with figures 3. The 
regions for which a zero of a occurs for $ ^ are shown with a gray shadow, while the blank regions are "forbidden" areas 
corresponding to zeroes of 3>. The values of A; ' , A^ ' in the gray areas comply with conditions (37) and so satisfy conditions 
(31a), (31b) and (31c). The level curve log 10 (j/) = 2.4 has been qualitatively sketched. This curve is extremely close (~ 10 -8 ) 
to the diagonal line corresponding to |Aj ' = A ; , marked by the letter "A" (the letter "B" marks the line A, ' = — A^ . 
By following the analogy with Perturbation Theory (see section 8), the line "A" can be associated with "adiabatic" initial 
conditions, while "isocurvature" initial conditions would be characterized by the horizontal axis Aj ' = 0. However, the values 
of y = y D for the latter condition are much smaller than y — 10 2 ' 4 for all \a[ ' | / (see figure 3b). Hence, initial conditions 
with A 2 - s ' / and leading to the right decoupling temperature (close to the level curve log 10 (y) = 2.4) cannot be termed 
"isocurvature", but rather "quasi-adiabatic" initial conditions. 



FIG. 5. The decoupling temperature T — T D = 4 x 10 3 K for y = y D . This figure displays the implicit plot of equation 
(43) in terms of A s , Aj and log 10 (y). The values of y for this temperature value are clearly shown to be very close to 
y = y D = 10 2 - 4 . 



FIG. 6. The equation a — for y = y D . Implicit plot of the solution to the equation a = in terms of log 10 (|A^ |) and 
log 10 (|A^ 'l), denoted as "log(Ds)" and "log(Dr)". This plot illustrates the constraints on the orders of magnitude of these 
quantities under the condition that a zero of a occurs for y — y D « 10 2 ' 4 . For values of |A ; ' | w 10 -4 , compatible with observed 
anisotropy of the CBR, we have the following estimated value: |A* S '| w 10 -8 . This can be a constraint on entropy fluctuations 
in primordial quasi adiabatic perturbations. 



FIG. 7. Comparison between r, the timescales of Thomson (t 7 ) and Compton (t c ) scatterings and t H . The figure 
depicts the ratios \og w (r/t H ), log 10 (t c /t H ) and log 10 (t 7 /t H ) along the range < log 10 (y) < 3 for A' m) = (3/4)A* r) — A* s) 
with Af = 10- 4 and 1(T 8 < A< s) < 10- 4 . The curves that branch out correspond to log 10 (r/t H ) for the displayed values of 
A^ s '. The thick curves are \og 10 (t c /t H ) (above) and log 10 (t 7 /t H ) (below) for the same values of A' s ' and A^ r '. Notice that r is 
much more sensitive to changes in A^ s ' than t 7 and t c . It is clear from the figure that r w t 7 for A' s ' w 1CP 8 for temperatures 
near the decoupling temperature, obtained from the condition t H — t 7 and closely matching T D « 4 x 10 3 K. (marked by 
log 10 (y) w 2.4). The figure also reveals that r and t c are of the same order of magnitude for higher temperatures closer to 
t = U. This is consistent with the fact that Compton scattering is the dominant radiative process in this temperature range. 
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Delta [m ] 



2.44 



2.42 



log(y)2.4 



2.38 



2.36 




